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Abstract 

In 1994 we showed that very large classes of systems of nonlinear PDEs 
have solutions which can be assimilated with usual measurable func- 
tions on the Euclidean domains of definition of the respective equa- 
tions. Recently, the regularity of such solutions has significantly been 
improved by showing that they can in fact be assimilated with Haus- 
dorff continuous functions. The method of solution of PDEs is based 
on the Dedekind order completion of spaces of smooth functions which 
are defined on the domains of the given equations. In this way, the 
method does not use functional analytic approaches, or any of the cus- 
tomary distributions, hyperfunctions, or other generalized functions. 



Type independent existence and regularity results for large 
classes of systems of nonlinear PDEs 

Ten years ago, in [4], the following significant threefold breakthrough 
was obtained with respect to solving large classes of nonlinear PDEs, 
see MR 95k:35002. Namely : 
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a) arbitrary nonlinear PDEs of the form 



(1) F(x, U(x), DPU(x), ...) = f(x), xett 

with F jointly continuous in all it arguments, / in a class of measur- 
able functions, Q C R n arbitrary open, p G N™, with \p\ < m, for 
m G N arbitrary given, and the unknown function U : ft — ► R, were 
proven to have 

b) solutions U which can be assimilated with usual measurable func- 
tions on Q, and 

c) the solution method was based on the Dedekind order completion 
of suitable spaces of smooth functions on Q. 

In fact, the conditions at a) can further be relaxed by assuming that 
F may admit certain discontinuities, namely, that it is continuous 
only on (Q \ E) x R m *, where S is a closed, nowhere dense subset 
of Q, while m* is the number of arguments in F minus n. This re- 
laxation on the continuity of F may be significant since such subsets 
of discontinuity £ can have arbitrary large positive Lebesgue measure. 

The method of order completion and the results on the existence and 
regularity of solutions can easily be extended to systems of nonlin- 
ear PDEs of the above form (1). Furthermore, initial and/or bound- 
ary value problems can be dealt with easily by this order completion 
method. 

In this way, the solutions of the unprecedented large class of nonlinear 
PDEs in (1) can be obtained without the use of any sort of distribu- 
tion, hyperfunctions, generalized functions, or of methods of functional 
analysis. Moreover, one obtains a general, blanket regularity, given by 
the fact that the solutions constructed can be assimilated with usual 
measurable functions on the corresponding domains Q in Euclidean 
spaces. 

Recently, in collaboration with R. Anguelov, see [1], a further signifi- 
cant improvement of the above mentioned 1994 results was obtained. 
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Namely, this time we can further improve the regularity properties of 
the solutions by proving that they always belong to the significantly 
smaller class of Hausdorff continuous functions on the open domains 
Q, see Appendix for a short account on Hausdorff continuous func- 
tions. 

It should be noted that the results in [4] on existence of solutions do 
for the first time in the literature manage fully to overcome the cele- 
brated 1957 Hans Lewy impossibility see [5], and in fact do so with a 
large nonlinear margin. 

Also, the existence results in [4], and thus their mentioned recent im- 
provement with resepct to the regularity of solutions, when solving 
large classes of nonlinear PDEs, supersede to a good extent the earlier 
similar ones obtained through the algebraic nonlinear theory of gener- 
alized functions introduced by the author in the 1960s, and developed 
since then alone or in collaborations, see 46F30 in the AMS Subject 
Classification at www.ams.org/index/msc/46Fxx.html, as well as [6], 
[7], [3, p. 7] and the literature cited there, or MR 89g:35001, MR 
92d:46098, Zbl.Math.717 35001, Bull.AMS, Jan. 1989, Vol. 20, No. 
1, 96-101. 

To further facilitate the understanding of the above mentioned results, 
it may be useful to point to the following. In his latest 2004 edition 
of his Springer Universitext book "Lectures on PDEs", see [2], V I 
Arnold starts on page 1 with the statement : 

"In contrast to ordinary differential equations, there is no unified the- 
ory of partial differential equations. Some equations have their own 
theories, while others have no theory at all. The reason for this com- 
plexity is a more complicated geometry ..." (italics added) 

However, as the above mentioned results show, since [4], there is an 
existence and regularity of solutions theory for the large class of sys- 
tems of nonlinear PDEs of the form in (1). Moreover, recently, the 
regularity result has been significantly improved by proving the ex- 
istence of Hausdorff continuous solutions for such general nonlinear 
systems of PDEs. 
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Appendix. Hausdorff continuous functions 

We shall deal with functions whose values can be usual or extended 
real numbers, that is, elements in 

R = RU {-oo, +00} 

Moreover, we shall allow the values of the functions to be not only 
numbers in R, but also closed intervals of such numbers, namely 

[a, b] C R, a, b G R, a < b 

It turns out to be quite surprising how much more appropriate such a 
framework is when one deals with large classes of nonsmooth functions 
in what is usually called Real Analysis. 

Indeed, by considering such interval valued functions one obtains a 
systematic and effective way to study and deal with a large variety of 
nonsmooth functions. Furthermore, one can gain important insights 
into the properties of such nonsmooth functions, properties which in 
fact are not available in the usual approach. 

It appears therefore that with the emergence in the second half of 
the 19-th century of a rigorous approach to Analysis, and specifically, 
with the Dirichlet definition of a function as having values given only 
and only by one single number, a certain undesired limitation was 
imposed in an unintended manner, especially what the study of non- 
smooth functions is concerned. 

Somewhat later, towards the end of the 19-th century, Baire brought 
in the concepts of lower and upper semi- continuous functions, when 
dealing with nonsmooth real valued functions. And in effect, he asso- 
ciated with each real valued function /, two other real, or extended 
real valued functions /(/) and S(f), with /(/) < f < S(f), which 
proved to be particularly helpful. However, following the prevailing 
mentality, each of these three functions were considered as being single 
valued. 
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As it turns out, however, by considering interval valued functions, such 
as for instance F(f) = [1(f), S(f)], one can significantly improve on 
the understanding and handling of nonsmooth functions. 
The study of interval valued functions can, among others, show that 
the particular case of functions which have values given by one single 
number is appropriate for continuous functions only. On the other 
hand, nonsmooth functions are much better described by suitably as- 
sociated interval valued functions. 

Indeed, in the case of functions / which are not continuous, a much 
better description can be obtained by considering them given by a pair 
of usual point valued functions, namely / = [/,/], thus leading to 
interval valued functions. And then, a natural class which replaces, 
and also extends, the usual point valued continuous functions is that 
of Haus dor ff- continuous interval valued functions. The distinctive and 
essential feature of these Hausdorff-continuous functions / = [ /, / ] 
is a condition of minimality with respect to the gap between / and /, 
with the further requirement that / be lower semi- continuous, and / 
be upper semi-continuous. 

In retrospect, it is surprising to see how near Baire came to such 
a treatment of nonsmooth functions, what deep results he obtained, 
and its correspondent for lower semi-continuous functions, and yet 
followed the prevailing trend which considered functions as having to 
have point, and not interval values. 

A good measure of the naturaleness of interval valued functions can 
be seen in the results related to the Dedekind order completion of var- 
ious spaces of continuous functions. And it is precisely such recently 
obtained results which allow for the mentioned significantly increased 
regularity properties of solutions of PDEs. 

These Dedekind order completions prove to be subspaces of Hausdorff- 
conti- 
nuous, thus interval valued functions. By the way, the space of Hausdorff- 
continuous functions itself is order complete. 

Such results extend easily to functions defined on large classes of topo- 
logical spaces. 

A further indication of the natural role interval values play in the study 
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of nonsmooth functions can be found in the Differential and Integral 
Calculus being presently developed for functions with such values. 

It will be useful to start by introducing a few notations. Let 

(A.l) IR = { [a, a] \ a, a G R = RU{— oo,+oo}, a < a} 

be the set of all finite or infinite closed intervals. 
The functions which we consider can be defined on arbitrary topo- 
logical spaces Q. For the purposes of the nonlinear PDEs studied in 
this book, however, it will be sufficient to assume that Q C R n are 
arbitrary open subsets. 

Let us now consider the set of interval valued functions 



(A2) A(n) = {f-.n — ► IR } 

By identifying the point a G R with the degenerate interval [a, a] G 
IR, we consider R as a subset of IR. In this way A(Q) will contain 
the set of functions with extended real values, namely 

(A3) A(tt) = {f:Sl — > R } C A(n) 

We define a partial order < on IR by 

(A4) [a, a] < [b, b] <^ a < b, a < b 

Now on A(f2) we define the partial order induced by (A2.1.4) in the 
usual point- wise way, namely, for /, g G A(f2), we have 

(A5) f<g ^ f(x)<g(x), x£ll 

Clearly, when restricted to A(Q), the above partial order on A(Q) 
reduces to the usual one among point valued functions. 

Given an interval a = [a, a] G IR, we denote 
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a — a if a, a finite 

oo if a = oo and a finite, or 
a = — oo and a finite, or 
a = — oo and a = oo 

if a — a — ±00 
which is called the width of the interval a. Also, we denote by 

(A7) \a\ = max{ \a\, \a\ } 

the modulus of the interval a = [a, a] G IR. 

In this way 

(A8) A{Q) = { f e A(0) I w(f(x)) = 0, x e Vl } C A(fi) 
Let / G A(Q). For every rr G fi, the value of / is an interval, namely 
f(x) = [ l(x), J(x) ], with l(x), J(x) G R, l(x) < J(x) 
Hence, every function / G A(Q) can be written in the form 

(A9) / = [/, / ], with /, 7e^(n), / < / < 7 

and 

(A10) / G A(Sl) ^ f = f = f 

In the particular case of functions in A(fl), that is, with extended real, 
but point, and not nondegenerate interval values, a number of results 
in the sequel were obtained by Baire. 

Most of the more general results concerning functions in A(fi), that 
is, with values finite or infinite closed intervals, have recently been 
developed by Anguelov. 

For x G Q, we denote by V x the set of all neighbourhoods V C Q of x. 
Let us consider the pair of mappings /, S : A(Q) — > A(Q), called 



(A6) w(a) = 
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lower and upper Baire operators, respectively, where for every func- 
tion / G A(f2), we define 

(All) I(f)(x) = sup VeVx inf {zef(y) \ y eV } 

(A12) S(f)(x) = mf V£Vx sup {zef(y) \ y eV } 

In Baire, these two operators were considered and studied in the par- 
ticular case of functions / G A(Vl). 

In view of the main interest here in this book in interval valued func- 
tions / G A(f2), it is useful to consider as well the following third 
mapping, namely, F : A(Q) — > A(f2), defined by 

(A13) F(f)(x) = [I(f)(x), S(f)(x)\, /GA(fi), xen, 

and called the graph completion operator. 

The lower and upper Baire operators / and S, and consequently, the 
graph completion operator F, applied to any interval valued function 
/ = [ f i f ] £ A(f2) can now be conveniently represented in terms of 
the functions / and /. Indeed, from (A2.1.11), (A2.1.12) it is easy to 
see that 

(A14) /(/) = /(/), S(f) = 

Hence F(f) can be written in the form 

(A15) F(f) = [1(f), S(f)} = [1(1), S(7)] 

Let us note that for every function / G A(Q) we have the relations, 
see (A5), (A9) 

(A16) /(/) = /(/) </</</< S(f) = S(f) 

and, thus, the inclusions 

(A17) f(x)CF(f)(x), xeQ 
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Furthermore, the lower Baire operator / : / — > /(/), the upper 
Baire operator S : f — > £>(/) and the graph completion operator 
F : f — > = [ /(/), £(/) ] are all monotone with respect to the 

order < in (A5) on A(Q), which means that for every two functions 
/, g e A(J7) we have 

(A18) /< 5 =► 1(f) < 1(g), S(f)<S(g), F(f) < F(g) 

The operator F is also monotone with respect to inclusion, namely 

(A19) f(x)Cg(x), xeQ =}► F(f)(x)CF(g)(x), xeQ 

With an immediate extension of Baire, one can also show that all three 
operators are idempotent, that is, for every / e A(f2), we have 

(A20) /(/(/)) = /(/), S(S(f)) = S(f), F(F(f)) = F(f) 
Definition Al 

A function / e A(Q) is called segment-continuous, or in short, s- 
continuous, if and only if 

(A21) F(f) = f 

□ 

In view of (A17), it is obvious that condition (A21) is equivalent with 

(A21*) F(f)(x)Cf(x), xeQ 

Furthermore, (A20), (A21) give 

(A22) F(f) is s-continuous for / e A(fl) 

Example Al 
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Let us illustrate the concept of s-continuity in the simplest case of 
functions with one variable and one single discontinuity. Thus, with 
Q = R, we take / : Q — > R, or in other words, / G A(£l), defined by 



/(*) = 



a if x < 
b if x = 
c if x > 



where a, b, c G R, a 7^ a Then / is not s-continuous. 

Let us now take / : Q, — > IR, that is, / G A(f2), defined by 



/(*) 



a if x < 

[ 6, c ] if x = 

ci if x > 



where a, b,c,d G R, a < d and b < c. Then f is s-continuous, if and 
only if 6 < a and d < c. Similarly, if a > d and b < c, then f is 
s-continuous, if and only if b < d and a < c. 

Consequently, returning to the first example above, it follows that f 
is s-continuous, if and only if a = b = c, that is, if and only if f 
is continuous, see (A2.1.32) below, for the general case of functions 

□ 



The fundamental concept is presented now in 



Definition A2 

A function / G A(f2) is called Haus dor ff- continuous, or in short, H- 
continuous, if and only if / is s-continuous, and in addition, for every 
s-continuous function g G A(f2), we have satisfied the minimality con- 



10 



dition on / : 

(A23) g{x) C f(x), xeVl 



9 = f 



We shall denote by H(Q) the set of all Hausdorff-continuous interval 
valued functions on Q. 



Example A2 

Let us again consider the second function in Example Al above. Then 
/ is H-continuous, if and only if a = b and c = d 

Let us give three further examples. 

First, let us define a : R — ► IR by 
- 1 if x < 

<!(.)•)= [ — 1, 1 ] if x — 

1 if x > 

and then we can define (3 : R 2 — ► IR by 

a(sm(l/(x 2 + y 2 ))) if (x, y) ± (0, 0) 

[-1,1] if (x, y) = (0, 0) 



P(x, y) 



It is easy to see that both a and (5 are H-continuous. 

The third example is a typical shock wave solution of the well known 
nonlinear PDE in Fluid Dynamics 

U t + UU X = 0, t > 0, x e R 

which corresponds to the initial value problem 
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1 if x < -1 



U{0, x) 



x if - 1< x < 



if x > 

Namely, with = [0, oo) x R, we have the solution U : Vt 
given by 



1 if < t < 1, x < t - 1 

x/(t-l) if < t < 1, t-l<x<0 

if 0<*<1, x>0 

1 if * > 1, x < (t- l)/2 

[ - 1, 1 ] if t>l, x= (t — l)/2 
if t > 1, a; > (t- l)/2 



IR 



C/(t, x) = 



Then [/ is H-continuous. 



Remark Al 

The minimality condition (A23) in the above definition of H-continuous 
functions proves to play a fundamental role. 

□ 

As for the significance of the regularity property of being Hausdorff 
continuous, here we an important similarity between usual continuous, 
and on the other hand, Hausdorff-continuous functions, on the other. 
Namely, both of them are determined uniquely if they are known on a 
dense subset of their domains of definition. 

This property comes in spite of the fact that Hausdorff-continuous 
functions can have discontinuities on sets of first Baire category, and 
such sets can have arbitrary large positive Lebesgue measure. 
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Indeed, we have 



Theorem Al 

Let f — [ f, f ], g = [ g, ~g ] e A(f2) be two H-continuous functions, 
and suppose given any dense subset D C fl Then with the partial 
order in (A5), we have 

a) l(x) < g(x), xeD =>- f < g on Q 

b) J(x) < g(x), xeD f < g on n 

c) f(x) < g(x), xeD =^ f < g on Q 
Also 

d) /(x) = ^(x), xefl ^ f = g on Q 
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